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ABSTRACT 

We have observed 37 bright, polarized radio sources with the Allen Telescope Array 
(ATA) to present a novel analysis of their Faraday rotation properties. Each source was 
observed during the commissioning phase with 2 to 4 100-MHz bands at frequencies 
ranging from 1 to 2 GHz. These observations demonstrate how the continuous frequency 
coverage of the ATA's log-periodic receiver can be applied to the study of Faraday 
rotation measures (RMs). We use KM synthesis to show that wide-bandwidth data can 
find multiple KM components toward a single source. Roughly a quarter of the sources 
studied have extra RM components with high confidence (brighter than ~ 40 mJy), 
when observing with a RM resolution of roughly 100 rad m~^. These extra components 
contribute 10%-70% of the total polarized flux. This is the first time multiple RM 
components have been identified in a large sample of point sources. For our observing 
configuration, these extra RM components bias the measurement of the peak RM by 10- 
15 rad m~^; more generally, the peak RM cannot be determined more precisely than the 
RM beam size. Comparing our 1-2 GHz RM spectra to VLBA polarimetric maps shows 
both techniques can identify complicated Faraday structures in the sources. However, 
the RM values and fractional polarization are generally smaller at lower frequencies than 
in the higher-frequency VLBA maps. With a few exceptions, the RMs from this work 
are consistent with that of earlier, narrow-bandwidth, all-sky surveys. This work also 
describes the polarimetry calibration procedure and that on-axis ATA observations of 
linear polarization can be calibrated to an accuracy of 0.2% of Stokes I. Future research 
directions include studying the time-dependent RM structure in Active Galactic Nuclei 
(AGNs) and enabling accurate, wide-area RM surveys to test models of Galactic and 
extragalactic magnetic fields. 
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Introduction 



Radio waves encode information not only about their origin, but about their entire path of 
propagation. One way this information is encoded is through Faraday rotation, the frequency- 
dependent rota t ion of the r adiation polarization angle caused by dispersion in a magnetized plasma 
(jFaradayl Il844l : iBurnI 119661 ) . Observationally, Faraday rotation is parameterized by the rotation 
measure (RM): 

RM = A0/A(A2), (1) 

where 9 is typically measured in radians and A in meters. The plasma dispersion law predicts RM 
induced by propagation along a line as: 



RM = 0.81 / UpB -dlradm" 



(2) 







where Ue is in cm ^ and B is in fiG, and d is the distance to the source ( BurnI 119661 ). Thus, 



measurements of Faraday rotation constrain physical conditions critical to understanding a wide 
variety of problems. 

Measurements of RM have expanded our knowledge of magnetic fields in our own Galaxy and 
in other ga laxies. Several efforts have been made to coin pile RM along lines of sight over large areas 
of the sky (jSimard-Normandinlll98ll : lBrown et al.ll2003l ). Finding pat terns in thes e RM values have 
been used to constrai n Galactic magnetic structure and turbulence (jHeilesI Il989l : iHan et al.l 120061 : 
Haverkorn et al.ll2008l ). Similar studies of oth er galaxies constrain models for the amplification of 
galactic magnetic fields (jGaensler et al.ll2005l ). Observations of individual AGNs and the massive 
black hole in the G alactic center have been used to constrain the geometry of their magnetic fields 



and depolarization (j Bower et al 



199S 



Zavala fc Tavloill2005l ). 



Recently, Taylor et al. (2009; hereafter "T09") , have expanded the number of sources with 
measured RM by nearly two orde rs of magnitude. T his was done by reanalyzing the 1.4 GHz 
NRAO VLA Sky Survey (NVSS; ICondon et al]ll998l l. producing RMs for 37,543 radio sources 
located throughout the sky north of declination -40°. Each source was observed in two bands, 
whi ch, a ssuming Equation [H can be used to estimate RM. The high density and large coverage of 
the lT09l sample make it very effective at statistically measuring the structure and strength of the 
magnetic fields. 



H oweve r, while the |T09| RMs are undoubtedly precise, it is not clear that they are accurate. 



BurnI (|l966l ) first noted that the Faraday rotation does not require the polarization angle to change 
as A^. When the emitting and Faraday-rotating media are mixed or multiple sources with dif- 
ferent RM are sp atially unresolved, one observes complicated changes in the polarization angle 
with wavelength (jGoldstein &: Reedl Il984l ). These kinds of changes, combined with the possibil- 
ity of mr ambiguities ii i mea suring polarization angle, make it difficult to measure RM robustly. 
Brent j ens &: de Bruynl ( 20051 ) noted that discrete sampling of the Stokes vector in A^-space con- 
strains different kinds of Faraday structures. They introduce the concept of "rotation measure 
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synthesis" , which takes advantage of the mathematical similarity between how aperture synthesis 
is done with multiple antennas and how RM is measured by sampling multiple wavelengths. The 
output from this technique is a "RM spectrum" showing the amou nt of polarized brightness as a 
function of its Faraday rotation (jHeald et al.ll2009l : iMao et al.ll2010l ). 



In parallel with this algorithmic development has been the technical development of wide- 
bandwidth radio receivers and digital signal processing. The Allen Telescope Array (ATA), a radio 
i nterferometer in n orthern California, is commissioning these and many other new technologies 
(jWelch et al.ll2009l ). One strength of the ATA design is its log-periodic receiv er, which gives it 
continuous access to frequencies from 0.5 to 1 GHz (IWilliams &: Bowed 120101 ). The array and 
receiver design are optimized for large surveys (jCroft et al.ll20ld ). which could be very powerful in 
the study of cosmic magnetism. 

The coincidence of these new algorithms and technologies for wide-bandwidth polarimetry 
inspired us to conduct a commissioning survey with the ATA. The first goal of this survey is to test 
the polarimetry capab ilities of the ATA telescope design. Second, comparing the RM measured by 
the ATA and by |T09| can test for biases associated with narrow bandwidth and coarse frequency 
resolution of the NVSS data. If the lT09l RM measurements are trustworthy, it will give confidence 
in their application to measuring large-scale RM structure. Finally, the wide bandwidth and fine 
frequency resolution of the ATA demonstrates the power of RM synthesis to study complicated 
RM spectra. The paper begins with a description of the observations in ^ The data reduction, 
including a detailed description of the ATA polarimetry calibration process, is described in f|3l The 
analysis of the RM spectra is given in ^ and implications of this work are given in ^ 



2. Observations 

The ATA is currently composed of 42 6.1-m dishes with an offset Gregorian optical design 
on an alt-az mount. The dishes are distributed to produce a Gaussian distribution of baselines 
with a maximum length of 300 m (1500A at 1.4 GHz). Typical observations will have a resolution 
of 150" at the zenith and a primary beam FWHM of 2?5 at 1.4 GHz. As shown in Fig. [H the 
receiver has a pyramid sh ape with a log-periodic pattern, giving it sensitivity from 0.5 to 10 GHz 



(jWilliams &: Bowed l2010l ). The entire bandpass is sent as an analog signal over ~100m of optical 
fiber into the on-site signal processing room. In the signal processing room, the signal is mixed, 
filtered, and digitized before being passed to the backends. 

For our experiment, the signal was processed simultaneously by two correlators tuned to dif- 
ferent frequencies. Each correlator takes 32 dual-linear polarization antenna inputs and produces 
1024-channel spectra covering 100 MHz. Each visibility cross product {xx, xy, yx, yy) was inte- 
grated for 10 seconds. All fields were observed with the correlators tuned to 1.43 and 2.01 GHz, but 
several fields were also observed with correlators tuned to 1.0 and 1.8 GHz. In both cas es, the focus 



was set to the higher frequency, which gives the array good sensitivity in both bands (jWelch et al 
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Fig. 1. — The log-periodic receiver of the ATA. This picture shows the open antenna/receiver 
configuration used during commissioning. In the current configuration, the c-clamp is removed, the 
receiver is enclosed from below by a shroud, and from above by a radio-transparent radome. 
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2009). These frequencies were chosen as a rough "minimum redundant spacing" in A -space (jHeiles 



20071 ). which is important for creating good RM spectra. The frequencies were also selected such 



that some portion of their bands were clear of radio frequency interference (RFI). 



With the goals of testing|T09| and demonstrating the power of RM synthesis, we soug ht a t arget 



list of a few tens of sources. We started with the catalog of polarized sources presented by lT09l from 
the 1.4 GHz NVSS data. We built a list of 46 sources by selecting sources with linearly-polarized 
brightness greater than 200 mJy beam~^ and declinations greater than —30°. This list excludes 
two sources in lT09l that were not real sources (sidelobes of Gas A). Filtering on polarized brightness 
produced a source list with a range of locations and RM that are relatively easy to detect. Note that 
the target selection criteria include our calibrator, 3G 286. We included 3G 286 as a target in our 
survey, which we refer to as J133108-I-303032, as a test of our calibration procedure. Galibration 
is derived from specific observations of 3G 286 that exclude J133108-I-303032. Observations of 
J133108-I-303032 are treated as any other target, so we can use errors seen in J133108-I-303032 to 
estimate calibration errors toward all targets. 

The survey was conducted over three separate days, as summarized in Table [TJ In total, 44 
fields were observed over 34.5 hours at two or four frequencies. The typical spatial resolution at the 
zenith is shown for each frequency. The schedule was designed to observe each target with two, 7- 
minute scans at different hour angles to improve image quality; the actual observations range from 
one to six scans. On day 1, 41 fields were observed, eight of which had data collection problems 
that led us to conservatively reject them. Another two of the fields (centered on J195904-I-393805, 
J200259-I-410150) were empty and are likely to be sidelobes of Gassiopeia A in the NVSS catalog. 
On day 2, we observed eight new fields at 1.43 and 2.01 GHz and 11 fields from day 1 at 1.0 and 
1.8 GHz. Of these, three and two fields at 1.0 and 1.8 GHz, respectively, were too short and had 
uv coverage too poor to be useful. On day 3, the target list was a subset of 14 fields from day 2, 
observed only at 1.0 and 1.8 GHz. Of this list, one field was rejected for having poor data quality. 
After these data quality checks, we had 37 fields with at least two 100-MHz bands, as shown in 
Table M 



Table 1. Summary of Observations 



Date Start (UT) Length (hrs) Resolution (arcsec) Fields 

1.0 GHz 1.43 GHz 1.8 GHz 2.01 GHz 



19-20 Sep. 2009 19 20 ... 250 x 120 ... 180 x 90 41 

27 Nov. 2009 15 9 360 x 170 250 x 120 200 x 100 180 x 90 19 

4 Dec. 2009 18.5 5.5 360 x 170 ... 200 x 100 ... 14 
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Table 2. Fields Observed in ATA Polarimetry Survey 



T09 Source Common Name Integration time (minutes) 

1.0 GHz 1.43 GHz 1.8 GHz 2.0 GHz 



J005558+682218 


3C 27 




20 




20 


J005734-012258 


3C 29 




20 




21 


J010850+131831 


3C 33 




20 




20 


J012644+331309 


3C 41 




28 




28 


J022248+861851 


30 61.1 




28 




28 


J030824+040639 


30 78 




28 




28 


J035232-071104 


30 94 




14 




18 


J052109+163822 


30 138 




28 




28 


J063633-204233 






21 




20 


J074948+555421 


40 56.16 


7 


13 


7 


13 


J084124+705341 


40 71.07 


20 


34 


40 


34 


J094752+072517 


30 227 




18 




17 


J 104244+ 120331 


30 245 


7 


11 


7 


11 


Jl 13007-144927 


PKSJ1130-1449 




10 




10 


J122906+020305 


30 273 


7 


11 


14 


10 


J123039+121758 


M87 (jet) 


21 


7 


36 


7 


J123522+212018 






14 




14 


J12561 1-054720 


30 279 


14 




14 


11 


J133108+303032 


30 286 


21 


7 


31 


7 


J153150+240243 


30 321 


19 


14 


35 


14 


J160231+015748 


30 327 




22 




24 


J160939+655652 


30 330 




22 




14 


J162803+274136 


30 341 




7 




7 


J164258+394837 


30 345 




7 




7 


J165111+045919^ 




19 


20 


37 


20 


J165112+045917'' 






20 




20 


J172025-005852 


30 353 


7 


14 


42 


14 


J184226+794517 


30 390.3 


19 


7 


38 


7 


J192451-291431 


PKS1921-293 


20 


20 


20 


17 


J194114-152431 






21 




20 


J201713+334546 


40 33.50 




20 




20 


J211636-205551 


PKS J2116-2055 


31 


13 


37 


13 
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Table 2 — Continued 



T09 Source 


Common Name 


Integration time (minutes) 








1.0 GHz 1.43 GHz 1.8 GHz 


2.0 GHz 


J212344+250410^ 


3C 433 


20 


20 


J212345+250448^ 


3C 433 


14 


14 


J222547-045701 


3C 446 


13 


13 


J225357+160853 


3C 454.3 


20 


20 


J231956-272713 


PKS 2317-27 


14 


13 



^J165111+045919 and J165112+045917 are unresolved by ATA 
''J212344+250410 and J212345+250448 are unresolved by ATA 
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Data Reduction 



3.1. Flagging 



The visibility data in each band were edited to remove RFI and poorly-performing antennas. 
Since the array was being commissioned during this survey, some receivers were being refurbished 
and were not avai lable. System health tests revealed a few antennas with large xy beam squint 
([Heiles et al.ll200ll ) that we flagged entirely. Initial flagging removed 100 channels from each edge of 
the 1024-channel band, where sensitivity is low; this removes roughl y 20% of usable channels. Next, 
the ATA automated flagging routine, RAPID (iKeating et al.ll2009l ). was run on the data. RAPID 
iteratively removed antennas, times, and frequency channels with large-amplitude deviations from 
their means. Daytime observations are sometimes affected by flux from the sun falling in distant 
sidelobes; for simplicity, we manually removed all baselines with uv distance less than 200A. Finally, 
a round of r nanual flagging was done by visually inspecting the data as (time, frequency) plots 
(described in lWilliams &: Bowed l20ld ). The number of antennas remaining after flagging changed 
with frequency and time, but typically data from about 21 antennas were used in the final analysis. 
The result of all flagging reduced final sensitivity by roughly 50%. 

Table [2] shows the total integration time for each field and frequency after flagging. The total 
time on a field ranges from 7 to 42 minutes, with a typical duration of 18 minutes. This produced 
15.1 hours of dual-polarization data on the target fields. The combined useful integration time for 
all sources and calibrators is about 18 hours, which is about 50% of the total scheduled time. 



3.2. Calibration 

The correlator outputs a complex visibility for each polarization, frequency, time, and baseline. 
The observed visibilities change with the distribution of sources in the field and the response of the 
telescope, both of which can change with polarization, frequency, time, and baseline. A general 
expression of how an astrophysical signal is observed as visibilities is given by the "measurement 



equation" formalism (iHamaker et al.lll996l ). The measurement equation treats the emitted signal 
as a vector that is modified by a series of matrix operations before being detected as visibilities. 
The equation can be written as v = Je, where e is the emitted signal as a Stokes vector (I, Q, U, 
V), J is the Jones matrix that describes the transformation of that vector during propagation and 
detection, and f is a vector of the observed complex visibilities {vxxi Vxy^ Vyx, '^yy)- 

The measurement equation formalism is useful because we can use the tools of linear algebra 
to simplify effects of signal propagation and detection. In the case of the dual-linear feeds used by 
the ATA, the Jones matrix can be written as an ordinary product of several matrices: 



J = G*D*C*P*S 



(3) 
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where G represents the antenna gain, D is the "leakage" of the two feed polarization into one 
another, C is the rotation of the feed relative to the mount, P is the rotation of the feed relative to 
the sky (the parallactic angle), and S transforms the Stokes vector into the xy coordinate system. 
The latter three of these ter ms (C, P, S) are ge ometric transformations of the electromagnetic 



wave and have been derived (IHamaker et al.lll996l ). The feed is oriented such that the horizontal 
polarization is x and the vertical polarization is y. By rotating the Stokes vector into the antenna 
coordinate system (applying C and P), we can see how the rotated Stokes vector appears as 
visibilities modified by the antenna gains a nd leakages . Exp anding these equations and dropping 
second-order terms produces the following (jSault et al.lll996l ): 



"xy 



"yx 



\ vyy ) 



( 



\ 



9x 



d2y) 



9xy 

-gyx{dly - dl^) 

9yy 



9xx 




~9yy 





9xy 
9yx 







'^9xy 
^9yx 





\ / Irot \ 
Qrot 
Urot 
J \ Vrot J 



(4) 



where gxx,xy,yx,yy are the products of the gains of any two antenna polarizations, and dix,iy,--- are 
the complex leakages of y into x for antenna 1, of x into y for antenna 1, etc.. These are the 
parameters solved during calibration. 

In general, the gains and leakages can change with time, frequenc}!!], antenna, and polarization. 
The quasar 3C 286 is an excellent calibrator for complex gains and leakages since it is bright (12- 
17 Jy from 1-2 GHz), isolated, unresolved by the ATA, strongly linearly-polari zed 10%), and 
widely used as a flux standard. Furthermore, 3C 286 has RM « rad m~^ (jRudnick &: Jones 



19831 ). which makes it very easy to identify frequency-dependent calibration problems. As a result, 
each observing period was scattered with repeated, on-axis observations of 3C 286. On day 1, 
3C 286 was not up during the entire observation, so 3C 138 was used to extend the time- and 
antenna-dependent gain calibration in time; 3C 286 was always used exclusively for leakage and 
frequency-dependent gain calibration. The f requency dependence of total intensity of 3C 286 is 
based on the flux scale of iBaars et al.l (119771 ) as parameterized by lOtt et al.l (|l994l ) . The linear 
polarization of 3C 286 is assumed to follow the VLA modej^. 



The data were calibrated using MIRIAD (polarimetry calibration described in ISault et al 



1991 



). To account for the frequency dependence of the leakages, the data were split into frequency 
segments of 5 channels; most analysis is done with these ~0.5 MHz segments. Next, frequency- 
dependent gain calibration was done with MFCAL with a calibration interval of 1 hour, the typical 
time scale of phase variations for the ATA. Finally, frequency-independent gains and leakages were 



^Traditionally, the frequency dependence of the gain is treated separately as the "bandpass"; here we include 
frequency dependence to define a more general gain. 

^See |http : / /www . via . nrao . edu/ astr o/calib/manual/polcal . html [ R. Perley & N. Killeen, private communica- 
tion 
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calculated with GPCAL with a calibration interval of 1 hour. Equation [3] has four equations and 
eight unknowns, so a single measurement of the calibrator is not enough to solve both gains and 
leakages uniquely. One can use the change in parallactic angle to effectively provide multiple sets 
of {vxx, Vxy, Vyx, Vyy) aud (irot, Qrot ) Urot, ~^rot) for which there is one solution to the (time- 
dependent) gains and (constant) leakages. For our observations, GPCAL converged to a solution 
for an observation of 3C 286 over a range of parallactic angle greater than a few tens of degrees. 
Due to a scheduling error, observations of 3C 286 on day 3 did not have enough parallactic angle 
coverage to solve for the leakages (A0 ~ 10°). The day 3 gain and leakage solutions converged 
after more aggressively cutting antennas with known large (> 20%) leakage and using the leakages 
measured on day 2 as an initial guess for day 3. 

The final gain solution had an amplitude that changed by less than 10% and gain phase scatter 
of less than 10° within each day and band. After calibrating 3C 286, the gain and leakage solutions 
for each band, day, and frequency segment were applied to their corresponding target field data. 
The time-dependent gain solutions were applied to the target fields by linearly interpolating in 
time. 

The calibration and visualization of the data taught us several things. First, a gain amplitude 
calibration error is an important source of error in Stokes Q and U. This is the well-known challenge 
of measuring linear polarization with linear feeds (i.e., Stokes Q is formed by differencing two large 
numbers, the xx and yy visibilities). This kind of error limits us to using data within 1 hour of a 
gain calibration solution. Second, Stokes I can also mix with Q and U if one doesn't solve for the xy 
phase difference for all antennas. In this work, we know the Stokes parameters for 3C 286 and can 
solve for the absolute xy phase, which we assume does not change within a day. Third, not solving 
for the absolute leakages of all antennas will cause Stokes Q and U to mix, potentially biasing the 
measurement of RM. Fortunately, all of these parameters can be solved for with observations over 
a range of parallactic angle of a strongly linearly-polarized source like 3C 286. 

Figure [2] plots the x-polarization leakages for different bands, based on observations of 3C 286. 
The typical magnitude of the leakages is less than 10%, with a few antennas ranging up to 20%. 
The leakage changes in frequency by about 3% per 10 MHz with a tendency to make loops in (real, 
imaginary) space. Comparing repeated observations of 3C 286 shows that the leakages change by 
less than 5% over a time scale of a month. Observations with 3C 286 located off-axis find false 
Stokes (Q, U) signals typically less than 1% of Stokes I within a third of the half-power point of 
the primary beam (0?3 at 2.01 GHz). Applying calibration derived on-axis is likely reliable for the 
offsets up to 0?1 considered here. However, since the errors for off-axis sources are not explicitly 
measured, results for these sources should be treated caut i ously. The behayior o f leakage in the 



ATA is discussed in more detail elsewhere (ILaw et al.ll2010l ; Ide Villiers et al.ll2010l ). 
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Fig. 2. — Sample of the frequency-dependent polarization leakage of several antennas. Each dot 
shows the complex leakage in a 5 MHz segment of the band; later analysis uses a resolution of 0.5 
MHz. The lowest-frequency, 5-MHz segment is labeled with the antenna number. For clarity, only 
half of the antennas are shown. The plots show different bands: 1.0 GHz (top left), 1.43 GHz (top 
right), 1.8 GHz (bottom left), and 2.01 GHz (bottom right). 
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3.3. Calibrated Data Quality 

As a first-order check of data quality we inspected tlie calibrated Stokes parameters of simple 
sources. Figure [3] shows three plots that demonstrate the calibration quality at 1.43 GHz toward 
J225357+160853, a bright, compact source seen over a wide range of parallactic angle. The cal- 
ibrated visibilities for this compact source in (real, imaginary) space have Gaussian distributions 
centered at the expected flux densities. Plots of the Stokes parameters as a function of parallactic 
angle indicate leakage calibration errors, since those kinds of changes move with the feed as it 
rotates on the sky. The plot shows a change in Stokes Q, U, and V with parallactic angle by at 
most 0.1 Jy (less than 1% of Stokes I). Errors of this scale should be expected in any 0.5 MHz 
segment; as described in ^3A\ the error analysis accounts for this by measuring the noise in Stokes 
V images. 

The ultimate test of the data quality was to measure the Stokes Q, U and RM of 3C 286 when 
treated as a calibrator and as a target. As described earlier, we observed 3C 286 for calibration, 
which required long observations over a range of parallactic angle. But we also observed 30 286 
briefly as a target, sometimes called J133108+303032 (see Table[2]). The general calibration solution 
is derived only from observations of 3C 286 as a calibrator, which is transferred to J 133108+303032 
and all other targets. This approach allows us to quantify two kinds of error. First, when treating 
3C 286 as a calibrator, the calibration method will ideally produce the fluxes in the model of 3C 
286. In practice, there are deviations that indicate some aspect of the ATA calibration model that 
is not correct (ionospheric Faraday rotation, unflagged RFI, second-order leakage corrections, etc.). 
Second, when treating 3C 286 as a target, any additional errors can be attributed to the application 
of calibration solutions to different times and observing conditions. Thus, measuring the RM of 
our calibrator tests the validity of our calibration model and estimates systematic error in the RM 
measured in the target fields. 

Table [3] shows the errors in observed polarization fraction and RM for 3C 286 when treated as 
a calibrator and as a target. The deviation of the calibrated data from the model tells us about the 
calibration applied to targets in the survey. Each row shows a combination of data from different 
days or bands that have calibration applied to them. For example, data from days 2 and 3 are 
merged in two cases because the same targets were observed on those days in those bands. In 
contrast, day 1 and 2 observations at 1.43 and 2.01 GHz are treated separately because different 
targets were observed on both of these days. The last two lines of Table [3] are the most relevant, 
since they show errors measured for 3C 286 for combinations of bands used to generate the final 
Stokes Q and U spectra used to generate the RM spectra. 

All error quantities are measured from observed Stokes Q and U values as a function of fre- 
quency. The change in the polarization angle as a function of is used to measure RM (in practice, 
RM synthesis is used; see §4.ip . The observed Stokes parameters are compared to the expected 
frequency-dependence of 3C 286 to measure fractional polarization error, Perr- The absolute po- 
larization error is the difference between the modeled and observed polarization vector averaged 
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IQI V 225357+160853 1,4077 GHz 



■a, ° 3^ 



Stokes Parameters vs Parallactic Angle for J225357+160853 




-30 -20 -10 
Parallactic angle (deg) 



Fig. 3. — Top: Visibilities in (real, imaginary) space for J225357+160853 at 1.4 GHz after applying 
the 3C 286 calibration solution. Visibilities were averaged over a 0.5 MHz segment and 20 minutes 
of integration time. Stokes I has larger scatter because there is an unpolarized, ~ 2 Jy source off 
axis. Bottom: The Stokes I, Q, U, and V of J225357+160853 plotted as a function of parallactic 
angle at two frequencies. Stokes parameters from typical 0.5 MHz segments in the 1.43 and 2.01 
GHz bands are shown. The mean and standard deviation in Stokes parameters are measured from 
the distribution of values over all baselines and channels in the segment. The plots show that there 
are no significant changes in the Stokes parameters with parallactic angle. This plot also shows 
that the daily variation in the ionosphere is not significant. The data points with parallactic angle 
of 50° corresponds to sunset at the ATA, which should have the largest RM contribution from the 
ionosphere, but shows little deviation from the other points. 
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Table 3. Calibration Error toward 3C 286 Treated as Calibrator and Target 



Day, Frequency 


3C 286 as 


Calibrator 


3C 286 


as Target 




\Pcrr 1 


RM 


|Pcrr| 


RM 


-, (GHz) 


(%I) 


(rad m~^) 


(%I) 


(rad m~^) 


Day 1, 1.43 


0.26 + 0.02 


0.05 + 0.33 


0.16 + 0.06 


-3.22 + 1.12 


Day 1, 2.01 


0.18 + 0.03 


-0.09 + 2.30 


0.15 + 0.26 


92.5 + 29.12=^ 


Day 2, 1.43 


0.30 + 0.05 


-0.82 + 1.10 






Day 2, 2.01 


0.21 + 0.08 


1.32 + 5.17 






Day 2+3, 1.0 


0.27 + 0.02 


-1.18 + 0.27 


0.65 + 0.10 


-2.79+1.36 


Day 2+3, 1.8 


0.07 + 0.01 


1.73 + 0.24 


0.42 + 0.08 


0.75 + 3.79 


Day 1, 1.43-2.01 


0.20 + 0.02 


0.02 + 0.08 


0.15 + 0.08 


-0.76 + 0.34 


All"^ 


0.08 + 0.01 


0.02 + 0.02 


0.20 + 0.05 


0.56 + 0.08 



^Flagging reduced the bandwidth and RM constraint significantly for all 
2.01 GHz observations on Day 1. 

'^Day 1 data at 1.43 and 2.01 GHz is combined with day 2 and day 3 data 
at 1.0 and 1.8 GHz in this measurement. There is no significant difference in 
the result if we use day 2 data at 1.43 and 2.01 GHz instead. 
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over all frequencies. The error in this quantity is the standard deviation in the error vector for all 
channels in the band. This analysis shows that the error in the polarization vector expected toward 
targets is roughly 0.2% of Stokes I, equivalent to an absolute position angle error of ~ 1°. The 
RM error toward 3C 286, when treated as a calibrator, is generally consistent with the measured 
noise. The RM measured toward 3C 286, when treated as a target, has a slightly larger deviation 
from RM = 0, suggesting that systematic effects are present at the ~ 2a level. The RM measured 
between two and four bands (relevant for all target fields) is more reliable, but shows that the RM 
of the target fields are not more accurate than ~ 0.6 rad m~^. A second polarization calibrator, 
3C 138 (J052109+163822), was also calibrated as a ordinary target field and has a RM consistent 
with its known value (shown in 



The accuracy of the RM measurement is cons istent with the variation expected from the 



ionosphere, which is t ypically less than 1 rad m ^ ( MacMillan &: Mausll2005l : iBilitza &: Reinisch 



20081 : iMao et al.ll2010l ). Ionospheric Faraday rotation may exceed this limit at sunrise and sunset. 



which spans about three hours of observing on day 1 and 1 hour of observing on day 2. Figure 
[3] shows how the Stokes parameters change for one source during sunset, when ionospheric RM is 
expected to be largest. The constancy of Stokes Q and U in that case limits RM variation due to 
the ionosphere to less than 2 rad m~^. 



3.4. Imaging 



After applying calibration solutions to target fields, each frequency segment was imaged in 
all four Stokes parameters. Imaging parameters at each frequency were chosen to maximize image 
fidelity for measuring peak fiux density. In general, images were not forced to the same beam size, 
since most sources are unresolved by the ATA. As described in §4.2^ only one field (J123039+121758) 
was imaged with a fixed beam size because its RM spectra showed some beam-size dependence. 
Section [12] discusses two other sources, J172025005852 and J184226+794517, that are doubles only 
marginally resolved at 1.0 GHz. 

Visibilities for the target fields were Fourier transformed with a "robustness" of 0, which slightly 
downweights long baselines to improve image fidelity. Images at all wavelengths had pixel size of 10" 
and 1024 pixels on a side, giving a field of view of 2?8. The pixels overresolve the synthesized beam 
(ranging from 107" to 215" from 2.0 to 1.0 GHz) and the ima ges cover most o f the primary beam at 
all wavelengths (ranging from 1?8 to 3?6 from 2.0 to 1.0 GHz; iHuU et al.ll2010l ). helping make robust 
images. The region within 5 to 10' of the center of each image was "cleaned" by iteratively fitting 
and removing the brightest sources down to a threshold of 4.5 times the Stokes I image theoretical 
noise level. Inspection of the images showed no obvious signs of cleaning problems. All subsequent 
analysis is done on a restored image f ormed by convolving the best-fit sourc e distribution with an 
idealized, Gaussian synthesized beam ( Hogbom 1974 : Thompson et ahlboOll ). 



The radiometer equation predicts the Stokes I thermal noise for ATA should be a = SEFD/VBW*t, 
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where SEFD is the system equivalent flux density, BW is the bandwidth and t is the length of the 
observation. During these observations, we commonl y used 21 dishes with a typical system tern - 



perature of 65 K, which translates to SEFD= 450 Jy (jWelch et al.ll2009l : IWilhams k Bowed l20ld ). 
For a typical image of a 5-channel segment, BW = 0.5 MHz, t = 7 minutes, the ideal sensitivity is 
30 mJy beam^^. 



In practice, short ATA commissioning observations are dynamic range limited ( Croft et al 



20ld ). The dynamic range limit manifests as a featureless, noise-like signal for all Stokes parameters 
that scales with stokes I. Our 1.43 GHz Stokes I images have a maximum dynamic range of 400 and 
a typical value of about 150. ATA images of Stokes Q, U, V and polarized intensity tend to be less 
limited by dynamic range and have noise levels within a factor of 2-5 of theoretical expectations. 
However, fields with Stokes I brighter than 20 Jy have more noise in polarized intensity. This sets 
the ATA dynamic range limit for images of polarized intensity at 3000 relative to Stokes 10. 

Examples of the image quality in two fields are shown in Figure [H First we show an image of 
J133108-I-303032 to demonstrate the image quality for the calibrator field when it is treated as a 
target. As with all fields, the Stokes I image of J133108-I-303032 is dynamic range limited at a level 
of roughly 350; the noise in the polarized intensity image is about two times the ideal, so it does not 
seem to be dynamic range li mited. The second image shows field J123039-I-121758, which contains 



M87 (jRottmann et al.l 119961 ). As the brightest and most complex field in our survey, it represents 
the worst-case scenario for image quality. The Stokes I image shows negative bowls produced by 
fiagging visibilities with M?;-distance less than 200A. The polarized flux image shows that the M87 
core is depolarized and that the jet is the brightest polarized source in the fleld with a brightness of 
1.2 Jy beam~^. The image has a standard deviation of about 55 mJy beam~^, which corresponds 
to a polarization dynamic range of roughly 20 relative to the polarized peak and 2600 relative to 
the total intensity peak. So, while the polarized intensity image noise is large, it is largely caused 
by the extreme brightness of M87. 

The presence of calibration errors and a dynamic range limit for the ATA require all image 
error analysis to be based on the observed noise. Since these effects and thermal noise in Stokes 
Q, U, and V scale together, we can use Stokes V to estimate t he noise in Stokes Q an d U. True 



Stokes V emission is generally much less than 1% of Stokes I (jHoman &: War did 1 19991 ). which is 
too weak to be detected in our individual 0.5 MHz images. We inspected dirty images of Stokes 
V and found that most images were noise-like and suitable for estimating noise in Stokes Q and 
U. In some cases, calibration errors produced weak sources in Stokes V at a level less than 1%, 
as discussed in ^3.31 However, the uncleaned Stokes V images tend to spread this false emission 
and increase the standard deviation of pixel values in the image. For this reason, we set the error 



^Image dynamic range is defined here as the maximum of an image divided by its standard deviation measured 
off peak. 

''Comparing the polarized intensity noise to the peak polarized intensity is not as meaningful, since some sources 
are depolarized. 



13''33™ 33™ 31™ 30™ 29™ 13''32™30^ 33™0'' 30^ 31™0^ 30^ 30™0^ 30^* 39™D^ 30^ 

RA (J2000) RA (J3000) 

Fig. 4. — Images showing imaging quality for a typical and a challenging field averaged over the 
1.43 GHz band on day 1. Left: Gray scale shows linearly-polarized flux of field J133108+303032 
(a.k.a. 3C 286) with the gray scale bar at right given in units of Jy beam~^. Contours show the 
total intensity at levels of -0.25, 0.25, 0.75, and 2.25, 6.75 Jy beam~^. Right: Polarized and total 
intensity images of field J123039+121758. This field contains M87, the brightest and most complex 
source in this survey, to demonstrate the most challenging field studied in this survey. Contours 
show total intensity at levels of -1.3, 1.3, 3.9, 11.7, and 35.1 Jy beam~^. The core of M87 is very 
bright in Stokes I, but largely depolarized when averaging across the 1.43 GHz band. 
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for each image of Stokes Q and U to the standard deviation of pixel values in uncleaned Stokes V 
image. 



4. Analysis 



4.1. Rotation Measure Synthesis 

Traditionally, the RM was measured by fitting the standard, plasma dispersion law to 
measurements of 9 at multiple wavelengths. However, this technique does not account for ambiguity 
in the rotation of the polarization angle between sampling points in A^-space (known as the "nvr 
ambiguity" ). Also, in many cases the pol arization an gle does not follow the A^ law, giving misleading 
results from the angle- fitting technique (|Burnlll966l ). 



Brent i ens &: de Bruvnl (j2005l ) derived a more generally valid way of studying RM: rotation mea- 
sure synthesis. First, the technique treats the linear polarization as a complex vector (P = Q + iU) 
rather than a scalar {6). Now the polarization angle change is visualized as a rotation of the Stokes 
vector with frequency. Second, the concept of rotation measure, which is essentially an observed 
quantity, is changed to a more general "Faraday depth", (j). In this way, the polarization vector 
from each Faraday depth, F{(f ) \ can be integrated to produce the appa.rent c omplex polarization 
surface brightness (jBurnlll966l : iBower et al.lll999l : iBrentiens fc de Bruvnl l2005l ): 



(5) 



The similarity of Equation [5] to a Fourier transform makes it easy to invert and derive the 
distribution of polarized brightness in Faraday depth-space. Ind eed, the name "RM synthesi s" was 
chosen to show similarity of this technique to aperture synthesis (jBrentiens &: de Bruvnll2005l ) . The 
distribution of samples in A^-space is analogous to antennas in physical space; widely separated 
points constrain closely spaced RM and compact sources, respectively. The shortest wavelength 
determines the largest "Faraday thickness" detectable (A</> ~ "^/-^min)- ^ Fourier transform of the 
distribution of points in A'^-space shows the "rotation measure spread function" (RMSF), which is 
analogous to a "dirty beam" in aperture synthesis. 

The output from the RM synthesis algorithm is shown in Figure [5l This shows the analysis of 
all data toward 3C 286, when treated as a target (see Table [3]). The script takes Stokes (Q, U) and 
their error (see §3.4p as a function of frequency. The algorithm applies a Fourier transform to the 
complex Stokes vectors to generate a dirty RM spectrum. The RM spectrum is then interactively 
cleaned by fitting the RMS F to the dirty RM spectrum and subtracting the best fit component 
("RM-CLEAN" algorithm; llleald et allbood ). Since 3C 286 has only one RM component, the thin 
line in the bottom right plot shows a typical RMSF for an observation with 1-2 GHz data. 



Much as with image cleaning algorithms, RM-CLEAN has many parameters that control the 
depth and robustness of the final cleaned spectrum. We developed a two-step cleaning process to 
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avoid undercleaning or confusing RM sidelobes with noise. The initial RM cleaning threshold is 3 
times the mean of the noise of all input Q and U images. We then measure the standard deviation 
in the cleaned Q and U spectra for RM larger than 10 times the RM resolution; this is equivalent 
to measuring the noise in a cleaned image far from sources near the phase center. The final RM 
cleaning was stopped when the peak brightness was less than 3 times this measured noise level. 
The final RM spectrum is made by convolving the brightest 10 clean components with a Gaussian 
of width equal to the central part of the RMSF and adding in residuals left after cleaning. 



4.2. Polarized Source Detections 

Of the 37 fields, most contain only a bright, polarized source at the phase center. However, 
since the fields were chosen by a simple cut on 1.4 GHz polarized flux density, some fields are more 
complicated. Of the 10 complicated fields, five have a second, polarized source within 0?1 of the 
phase center, where the polarization calibration is believed to be reliable. 

The RM synthesis analysis was applied to both the central and off-axis sources. All complicated 
fields, including the coordinates of secondary sources, are listed in Table HI The RM synthesis 
analysis required a Stokes Q and U flux for each source at each frequency. For the sources at phase 



center, the fluxes were measured at the central pixel of each image, where the lT09l source is known 
to be. Two fields, J172025005852 and J184226+794517, have doubles resolved only at frequencies 
of 1.43 GHz and higher. This means that the 1.0-GHz polarized properties of the two sources in 
each field are mixed. Since these sources are observed with four bands, this ambiguity means that 
roughly a quarter of the total polarized flux may be shared between these double sources. 

Sources in fields J022248+861851, J063633-204233, and J123039+121758 were at least partially 
resolved. Since we measure fiuxes in Jy beam~^, the frequency dependence of the beam size could 
systematically affect flux measurements in these cases. For these three fields, we repeated RM 
synthesis analysis for these sources based on images with a fixed beam size at all frequencies. For 
sources in fields J022248+861851 and J063633-204233 there was no significant change to the RM 
synthesis results. In the analysis of field J123039+121758, there was significant difference, so all 
analysis is done on images with a beam size of 440" x 150" at a position angle of 74°. The sources 
in this field are discussed in detail in §5.4.41 



4.3. RM Spectra 

In total, we measured the Stokes Q and U spectra for 42 sources at a resolution of 0.5 MHz. 
Applying the RM-CLEAN algorithm to the Stokes Q and U spectra produced cleaned RM spectra 
for each source. The spectral resolution limits the RM spectra sensitivity beyond roughly ±90000 
rad m~'^. Cleaned RM spectra for all sources are shown in Figures [6] through 1461 
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133108+30: RM = 0.6 ± 0.1 rad/m^ 
Pol. Flux = 1319.36 mJy ; SNR = 261.74 




0.02 0.04 0.06 0.08 -4000-2000 2000 4000 

(Wavelength/m)^ RM (rad/m^) 



Fig. 5. — An example of the input to and output from the RM synthesis algorithm for all obser- 
vations of 3C 286, when treated as a target. Deviations of the RM clean model from the ideal 
(RM = 0) should be similar to the typical errors in RM spectra for all targets presented in this 
survey. Top right: Each point shows the measured Stokes Q and U brightness at one frequency. 
The clean model for the brightest 10 clean components is shown with a red line. Top left: The 
points show the measured polarization angle as a function of with a red line showing the clean 
model. Bottom left: The green, blue, and purple points show the measured polarized brightness. 
Stokes Q, and Stokes U, respectively, as a function of A^. The corresponding colored lines show the 
clean model. Bottom right: The RM spectrum from —5000 to 5000 rad m~^. The thin and thick 
lines show the dirty and clean RM spectra, respectively. 
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Table 4. Fields with Complexity 



Field Name(s) 


Morphology 


Second Source 




Stokes I 


Pol. Intensity 


RA, Dec (J2000) 


J010850+131831 


double 


double 


(1:08:55, 13:22:14) 


J022248+861851 


extended 


extended 




J063633-204233 


extended 


extended 




J123039+121758 


double 


extended 


(12:30:49, 12:23:23) 


J123522+212018 


double 


double 


(12:35:30, 21:20:48) 


J160231+015748 


double 


single 


(16:02:19, 1:58:25) 


J165111/2+045919/71 


single 


single 




J172025-005852 


double 


double 


(17:20:34, -0:58:43) 


J184226+794517 


double 


double 


(18:41:50, 79:47:28) 


J212344/5+250410/481 


single 


single 





Double resolved by VLA, unresolved by ATA 
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Fig. 6.— RM synthesis summary plots for 3C 27 (J005558+682218). {Top:) Plot showing the 
polarized flux (red stars) and polarization angle (blue points) as a function of A^. (Middle:) Plot 
showing the RM spectrum for —5000 < RM < 5000 rad m~^, derived from all data for each source. 
The dashed line shows 5 times the observed noise in the RM spectrum. (Bottom:) Plot showing 
the RM spectrum for -90000 < RM < 90000 rad m^^ derived from the 1.43 GHz data only. 
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Fig. 7.— RM synthesis summary plots for 3C 29 (J005734-012258), as in Figure El 



- 24 - 



010850+131831 



£! 



o 

CL 




0.020 

m 



0.025 0.030 0.035 0.040 

Lambda"2 {m^2) 



0.045 



0.050 




-50000 



50000 



RM Erad/m'^2) 



Fig. 8. — RM synthesis summary plots for 3C 33 (J010850+131831), as in Figure El 
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Fig. 9. — RM synthesis summary plots for the secondary source in the field with 3C 33 
(J010855+132214) as in Figure El 



- 26 - 



012644+331309 




-50000 



50000 



RM Erad/m'^2) 



Fig. 10.— RM synthesis summary plots for 3C 41 (J012644+331309), as in Figured 
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Fig. 11.— RM synthesis summary plots for 3C 61.1 (J022248+861851), as in Figured 
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Fig. 12.— RM synthesis summary plots for 3C 78 (J030824+040639), as in Figured 
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Fig. 13.— RM synthesis summary plots for 3C 94 (J035232-071104), as in Figured 
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Fig. 14.— RM synthesis summary plots for 3C 138 (J052109+163822), as in Figure El 
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Fig. 15. — RM synthesis summary plots for J063633-204233, as in Figure El 
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Fig. 16.— RM synthesis summary plots for 4C 56.16 (J074948+555421), as in Figure El 
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Fig. 17.— RM synthesis summary plots for 4C 71.07 (J084124+705341), as in Figure El 



- 34 - 



094752+072517 

90 F 1 1 1 




-50000 50000 

RM Erad/m'^2) 



Fig. 18.— RM synthesis summary plots for 3C 227 (J094752+072517), as in Figure El 
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Fig. 19.— RM synthesis summary plots for 3C 245 (J104244+120331), as in Figure El 
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Fig. 20.— RM synthesis summary plots for PKSJ1130-1449 (J113007-144927), as in Figure El 
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Fig. 21.— 



RM synthesis summary plots for 3C 273 (J122906+020305), as in Figure El 
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Fig. 22.— RM synthesis summary plots for the M87 jet (J123039+121758), as in Figure El 
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Fig. 23.— RM synthesis summary plots for the M87 core (J123049+122323), as in Figured 
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Fig. 24. — RM synthesis summary plots for J123522+212018, as in Figure [H 
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Fig. 25. — RM synthesis summary plots for the secondary source in the field of J123522+212018, 
as in Figure [H 
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Fig. 26.— RM synthesis summary plots for 3C 279 (J125611-054720), as in Figured 
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Fig. 27.— RM synthesis summary plots for 3C 321 (J153150+240243), as in Figure El 
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Fig. 28.— RM synthesis summary plots for 3C 327 (J160231+015748), as in Figure El 
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Fig. 29.— RM synthesis summary plots for 3C 330 (J160939+655652), as in Figure El 
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Fig. 30.— RM synthesis summary plots for 3C 341 (J162803+274136), as in Figure El 



-47- 



164258H-394S37 




0.020 0.025 0.030 0.035 0.040 0.045 0.050 

Lambda"2 {m^2) 




4000 



RM Erad/m'^2) 

Fig. 31.— RM synthesis summary plots for 3C 345 (J164258+394837), as in Figure El 
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Fig. 32. — RM synthesis summary plots for J165111+045919, as in Figure [H 
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Fig. 33. — RM synthesis summary plots for J165112+045917, as in Figure [H 
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Fig. 34.— RM synthesis summary plots for 3C 353 (J172025-005852), as in Figure El Note that 
J172025-005852 and J172034-005843 are unresolved at 1.0 GHz, so up to a quarter of the polarized 
flux may be shared between this and Figure [35l 
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Fig. 35. — RM synthesis summary plots for the secondary source in the field of 3C 353 (J172034- 
005843), as in Figured Note that J172025-005852 and J172034-005843 are unresolved at 1.0 GHz, 
so up to a quarter of the polarized flux may be shared between this and Figure [Ml 
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Fig. 36.— RM synthesis summary plots for 3C 390.3 (J184226+794517), as in Figure El Note 
that J184226+794517 and J184150+794728 are unresolved at 1.0 GHz, so up to a quarter of the 
polarized flux may be shared between this and Figure [371 
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Fig. 37. — RM synthesis summary plots for the secondary source in the field of J184226-I- 794517 
(J184150+794728), as in Figured Note that J184226+794517 and J184150+794728 are unresolved 
at 1.0 GHz, so up to a quarter of the polarized flux may be shared between this and Figure [36l 
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Fig. 38.— RM synthesis summary plots for PKS1921-293 (J192451-291431), as in Figured 
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Fig. 39. — RM synthesis summary plots for J194114-152431, as in Figure [H 
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Fig. 40.— RM synthesis summary plots for 4C 33.50 (J201713+334546), as in Figured 
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Fig. 41.— RM synthesis summary plots for PKS J2116-2055 (J211636-205551), as in Figure El 
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Fig. 42.— RM synthesis summary plots for 3C 433 (J212344+250410), as in Figure El 
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Fig. 43.— RM synthesis summary plots for 3C 433 (J212345+250448), as in Figure El 
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Fig. 44.— RM synthesis summary plots for 3C 446 (J222547-045701), as in Figured 
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Fig. 45.— RM synthesis summary plots for 3C 454.3 (J225357+160853), as in Figure El 
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Fig. 46.— RM synthesis summary plots for PKS 2317-27 (J231956-272713), as in Figure El 
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One goal of the observations was to find evidence for multiple RM components in these sources. 
To maximize our sensitivity to a wide range of RM components, we cleaned RM spectra in two 
ways. First, a high-resolution RM spectrum was made to measure components at small RM very 
precisely. Second, a high-sensitivity, low-resolution RM spectrum was created to search for large 
RM values. While there are more data and resolution when using the entire bandwidth range, 
spectral index and sampling effects make the RMSF less smooth; this is analogous to aperture 
synthesis with a sparse versus a dense array. 

For these reasons, we searched for components in RM spectra as shown in Figure [6] through 
H6l First, the full frequency range (either 1.43-2.01 GHz or 1.0-2.01 GHz) was used to generate 
RM spectra for —5000 < RM < 5000 rad m~^ with a resolution of roughly 51 or 141 rad m~^ (RM 



beam FWHM = 2\/3(A^ax - ^min); iBrentiens fc de BruvnIboOSi ). Second, the 1.43 GHz data 



which has the best quality and completeness, were used to generate RM spectra for —90000 < 
RM < 90000 rad m~^ with a resolution of roughly 688 rad Note that the largest Faraday 

thickness detectable is 150 rad m~^ for the high resolution spectra and 75 rad m~^ for the low 
resolution spectra. Thus, most of the RM spectra should observe unresolved RM components; 
the 1.0-2.01 GHz, high-resolution spectra are sensitive to mildly-resolved (~ 2 RM beam widths) 
Faraday structures. 



4.4. RM Components 

Table [5] gives the 5a threshold and RM resolution of each of the 42 RM spectra over the two 
RM ranges. As described in §4.H the noise in Stokes (Q, U) spectra was measured in cleaned 
RM spectra for RM larger than 10 times the RM resolution. The flux limit has units of mJy per 
synthesized beam per RM beam. The typical limits on the polarized brightness are roughly 40 mJy 
per 51 rad m"^ beam and 25 mJy per 141 rad m^^ beam for —5000 < RM < 5000 rad m^^. For 
-90000 < RM < 90000 rad m'^, the typical limit is 20 mJy per 688 rad m'^ beam. The last 
column shows the fraction polarization limit at 1.43 GHz over the —90000 < RM < 90000 rad m~^ 
range. 
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Table 5. Polarized Brightness Thresholds and Resolution of RM Spectra 

Source P±5ooo ^M^™ ^'igoooo ^M^Zoo^ /±90000 

(mJy/beam/RM beam) (rad m~^) (mJy/beam/RM beam) (rad m~^) (%) 



J005558+682218 


18 


141 


17 


688 


0.3 


J005734-012258 


22 


141 


21 


688 


0.5 


J010850+131831 


24 


141 


20 


688 


0.2 


J010855+132214'^ 


23 


141 


13 


688 


0.4 


J012644+331309 


21 


141 


13 


688 


0.3 


J022248+861851 


17 


141 


17 


688 


0.4 


J030824+040639 


19 


141 


14 


688 


0.2 


J035232-071104 


33 


141 


29 


690 


1.3 


J052109+163822 


12 


141 


11 


688 


0.1 


J063633-204233 


28 


141 


25 


688 


0.8 


J074948+555421 


59 


51 


136 


688 


10.5 


J084124+705341 


25 


52 


12 


688 


0.4 


J094752+072517 


30 


136 


36 


688 


0.5 


J104244+120331 


45 


51 


35 


688 


1.2 


J113007-144927 


31 


136 


29 


688 


0.6 


J122906+020305 


215 


51 


159 


688 


0.3 


J123039+121758 


117 


51 


122 


688 


2.0 


J123049+122323^ 


250 


51 


213 


688 


0.1 


J123522+212018 


22 


141 


17 


688 


1.1 


J123530+212048'' 


18 


141 


17 


688 


1.1 


J125611-054720 


87 


51 


101 


1523*^ 


1.1^ 


J133108+303032 


25 


51 


24 


688 


0.2 


J153150+240243 


18 


51 


19 


689 


0.6 


J160231+015748 


21 


141 


16 


688 


0.3 


J160939+655652 


22 


141 


17 


688 


0.3 


J162803+274136 


39 


141 


27 


688 


1.5 


J164258+394837 


35 


141 


23 


688 


0.3 


J165111+045919 


110 


51 


35 


688 


0.1 


J165112+045917 


120 


141 


33 


688 


0.1 


J172025-005852 


57 


51 


33 


688 


0.2 


J172034-005843^ 


55 


51 


44 


688 


0.1 
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Table 5 — Continued 



Source -P±5000 R-^iSOOO^ -P±90000 ^^'^±M)Oo¥ /±90000 

(mJy/beam/RM beam) (rad m~^) (mJy/beam/RM beam) (rad m~^) (%) 



J184226+794517 


76 


51 


25 


688 


0.4 


J184150+794728'' 


41 


51 


27 


688 


0.7 


J192451-291431 


77 


51 


34 


688 


0.4 


J194114-152431 


20 


141 


18 


688 


0.3 


J201713+334546 


9 


141 


11 


688 


0.4 


J211636-205551 


39 


51 


36 


688 


1.4 


J212344+250410 


17 


141 


19 


688 


0.2 


J212345+250448 


41 


141 


19 


688 


0.2 


J222547-045701 


39 


141 


37 


688 


0.6 


J225357+160853 


24 


141 


25 


688 


0.2 


J231956-272713 


38 


141 


42 


688 


2.0 



^Off-axis source. 

''No 1.43 GHz data available, so this limit measured at 1.8 GHz. 
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Tables [6] and [7] show the brightness, RM, and signal-to-noise ratio for components detected in 
the high- and low-resolution RM spectra, respectively. The component properties are measured by 
fitting a Gaussian to the cleaned RM spectra after subtracting the observed noise in quadrature. All 
RM peaks with a best-fit peak brightness greater than 5cr are shown. While we find several plausible 
candidates at this threshold level (see past experience suggests that some 5a candidates are 
not real. As such, the RM components in the tables should be considered candidates; most analysis 
presented here uses a stricter threshold of 7a. Varying parameters of the cleaning process (e.g., 
depth of initial clean, RM grid size, maximum RM) showed a roughly 20% variation in the brightness 
of components in a few of the most complicated RM spectra (noted in the tables). Aside from this 
caveat, all components shown in the tables are significantly detected regardless of the parameters 
used during cleaning. We also checked that RM spectra were robust to possible calibration errors 
in a few (Q, U) points (e.g., points with polarization angle > 45° in Figure [26j) . In all cases, the 
best-fit components were unchanged within their errors. 

In addition to being filtered based on a simple noise threshold, RM candidates were also fil- 
tered according to a RM dynamic range limit. The imaging and RM cleaning techniques described 
in ^3.41 and ^4.11 control for the effect of thermal noise and dynamic range limitations, but other 
systematic effects can create false RM components. Figure [5] shows how gain errors between the 
four bands create a wobble in the RM spectrum of 3 C 286, when treated as a target (absolute gain 



calibration accuracy w 3%; IWilliams &: Boweiil20ld ). This wobble requires two ~ 6a RM compo- 



nents with RM w —160 and 90 rad m ^ and polarized brightness of about 35 mJy. Furthermore, 
the spectral index of the sour ce can create false RM components by effectively modifying the RMSF 



(jBrentjens &: de Bruynll2005l ) 



The gain error and spectral index effects are similar in that they act as unexpected scalings 
of the polarized flux on frequency scales larger than the 100-MHz bandwidth. This kind of error 
will produce RM artifacts with polarized brightnesses that scale with the brightness of the source. 
This effectively creates a RM dynamic range limit for all sources with similar spectral indices and 
calibration errors. The RM spectrum of 3C 286 (treated as a target) allows us to estimate the RM 
dynamic range as the ratio of the true RM component to the largest false RM component. This 
estimate gives a dynamic range limit of 38 or, equivalently, 2.6% of the peak polarized brightness. 
This limit only removes the extra RM components seen in 3C 286. RM spectra with evidence 
for larger or more unusual spectral changes than 3C 286 (e.g.. Figure [23]) should be interpreted 
cautiously. 
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Table 6. Components Measured in High-resolution RM Spectra 





p 


RM 


SNR 






{ Vfici TTl 




J005558+682218 


305.6 


-99.0 ± 0.8 


84 


J005734-012258 


559.0 


0.5 ± 0.5 


123 


J010850+131831 


789.7 


-12.3 ± 0.4 


161 


J010855+132214- 


^32. 9 


-10.0 ± 1.9 


28 


J012644+331309 


238.5 


-68.0 ± 1.1 


56 


J022248+861851 


296.2 


-4.6 ± 0.8 


83 


J030824+040639 


153.9 


5.4 ± 1.6 


39 


J035232-071104 


191.6 


15.7 ± 2.0 


28 


J052109+163822 


622.4 


-0.7 ± 0.2 


249 


J063633-204233 


452.4 


45.5 ± 0.8 


78 


J074948+555421 


186.2 


1.5 ± 1.2 


15 




70.0 


54.4 ± 3.2 


5.9 


J084124+705341 


218.2 


-12.7 ± 0.6 


42 




43.1 


-161.1 ± 8.9 


8.3 




30.2 


-313.3 ± 11.4 


5.8 




32.3 


105.2 ± 39.2 


6.2 ^ 


J094752+072517' 


= 225.8 


-2.3 ± 1.6 


36 




47.8 


158.4 ± 12.2 


7.8 




34.2 


397.1 ± 81.0 


5.6 ^ 


J104244+120331 


165.0 


13.5 ± 1.7 


18 


Jl 13007-144927 


155.0 


42.5 ± 2.1 


24 




39.5 


-109.8 ± 5.6 


6.3 




32.8 


-242.3 ± 7.5 


5.2 


J122906+020305 


604.7 


-10.7 ± 2.3 


14 




277.3 


150.7 ± 4.1 


6.4 


J123039+121758 


306.4 


83.5 ± 2.3 


13 




190.7 


37.2 ± 3.1 


8.1 




137.0 


-118.1 ± 4.0 


5.8 




128.9 


-173.5 ± 3.9 


5.7 




209.5 


241.7 ± 2.6 


8.9 




136.5 


293.0 ± 4.3 


5.8 
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Table 6 — Continued 



Source 


P 


RM 


SNR 




fmJv) 


(rad m~^) 






129.5 


454.1 ±4.1 


5.5 




120.7 


507.0 ± 4.6 


5.1 


J123049+122323'' 


^=632. 2 


46.0 ± 1.4 


12 




560.8 


11.5 ± 2.2 


11 




295.4 


-66.2 ± 3.0 


5.9 




352.6 


101.2 ± 2.5 


7.1 


J123522+212018 


204.2 


-3.1 ± 1.4 


44 


J123530+212048^ 


95.4 


4.8 ± 2.6 


25 


J125611-054720 


283.5 


15.3 ± 1.1 


16 




165.9 


65.8 ± 1.7 


9.4 


J133108+303032 


1319.6 


0.6 ± 0.1 


261 


J153150+240243 


239.3 


11.4 ± 0.3 


63 




26.7 


77.1 ± 4.0 


7.1 




22.8 


-62.0 ± 4.0 


6.1 


J160231+015748 


399.9 


9.0 ± 0.7 


92 


J160939+655652 


242.3 


17.8 ± 1.2 


52 


J162803+274136 


281.6 


23.2 ± 1.8 


36 


J164258+394837 


368.2 


20.8 ± 1.2 


51 


J165111+045919'^ 


543.1 


-6.0 ± 1.4 


24 




134.0 


-127.0 ± 4.9 


6.1 




137.5 


107.5 ± 4.9 


6.2 


J165112+045917 


513.7 


-0.2 ± 2.5 


21 




125.9 


-152.2 ± 7.7 


5.2 


J172025-005852 


1065.4 


39.8 ± 0.3 


92 




158.6 


-68.7 ± 2.3 


13 




85.9 


134.4 ± 6.8 


7.5 




80.2 


196.5 ± 4.9 


7.0 


J172034-005843^ 


288.4 


107.5 ± 1.7 


25 




178.4 


-23.3 ± 2.9 


15 


J184226+794517 


354.6 


-8.4 ± 1.0 


23 




79.8 


-230.3 ± 4.0 


5.2 
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Table 6 — Continued 



Source 


P 


RM 


SNR 




(mJy) 


/ 1 9\ 

(rad m ) 




T1 841 iiD-U 704798^ 


ooo.u 


_l 4 c; -L n 

m u.u 


40 




Q1 8 


_-| 70 « _|_ 1 o 


1 1 




4Q Q 




u.u 




4^ 1 




^ 8 




zlQ 1 


_OOQ Q _L 


^ Q 




4fi 4 




5 fi 




4^ Q 




5 % 
0.0 








u.o 




51 1 


Q09 q _|_ 4 n 






4^ 


c^q 7 _|_ q 


5 9 




1 59 4 


1044- 1 Q 


Q 8 

iy.O 




1 3Q 


940 4 + 99 


Q 




1 95 


-49 4 + 90 


8 1 




Q4 7 


_q« _|_ 9 c; 


1 


71 QzLI 1 4—1 ^9zL'^l 


974 Q 


_7q 1 _|_ n q 


uu 




98fi 1 


000. 1 —1— U.t: 


144 


J211636-20555F 


198.3 


-22.6 ± 1.1 


24 




42.5 


185.4 ± 5.7 


5.3 




45.0 


-387.8 ± 5.6 


5.6 




46.7 


-236.9 ± 7.8 


5.8 




43.6 


-173.3 ± 10.0 


5.5 




41.6 


-104.0 ± 11.3 


5.2 


J212344+250410 


561.1 


-71.4 ± 0.4 


157 


J212345+250448 


440.1 


-71.7 ± 1.3 


52 


J222547-045701 


190.7 


-19.3 ± 3.0 


24 




105.9 


-149.9 ± 5.5 


13 


J225357+160853 


989.9 


-53.9 ± 0.3 


205 


J231956-272713 


165.7 


14.0 ± 2.9 


21 



-70- 



^Off-axis source. 

^Component is wider than the RM beam. 
Significance is underestimated. 

^RM component brightnesses have extra 
20% uncertainty introduced by RM cleaning 
algorithm. 
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Table 7. Components Measured in Low-resolution RM Spectra 



Source P 


RM 


SNR 


(ryi Tv I 


{ Vficl TTI 




J005558+682218 341.3 


-104.7 ± 3.1 


96 


J005734-012258 593.0 


0.9 ± 2.2 


137 


J010850+131831 861.5 


-14.6 ± 1.5 


205 


J010855+132214^126.1 


-3.3 ± 6.8 


45 


J012644+331309 253.9 


-66.6 ± 3. 


95 


18.8 ■ 


-1464.4 ± 55.4 


7.1 


J022248+861851 326.1 


-0.3 ± 3.2 


92 


J030824+040639 157.2 


9.5 ± 5.8 


52 


J035232-071104 191.0 


7.6 ± 9.5 


32 


J052109+163822 637.8 


-0.9 ± 1.0 


289 


J063633-204233 496.9 


51.4 ± 3.0 


99 


J074948+555421 217.4 


54.6 ± 41.5 


8.0 


J084124+705341 221.4 


-15.4 ± 3.4 


86 


J094752+072517 257.7 


21.4 ± 8.6 


35 


J104244+120331 156.2 


57.1 ± 14.0 


21 


J113007-144927 151.0 


22.2 ± 12.2 


25 


J122906+020305 778.0 


10.8 ± 12.7 


24 


J123039+121758 575.5 


75.4 ± 14.1 


23 


J123049+122323^218.8 


29.5 ± 59.0 


5.1 


J123522+212018 213.2 


-4.9 ± 4.8 


61 


J123530+212048^ 96.2 


-2.2 ± 11.4 


28 


J125611-054720 314.6 


232.7 ± 41.3 


15 


J133108+303032 1357.4 


-3.4 ± 1.1 


273 


J153150+240243 240.2 


-12.3 ± 4.9 


63 


J160231+015748 433.1 


-0.7 ± 2.2 


131 


J160939+655652 247.8 


17.6 ± 4.2 


71 


J162803+274136 297.9 


26.9 ± 5.5 


54 


J164258+394837 368.6 


20.4 ± 3.8 


78 


J165111+045919 435.2 


-29.5 ± 4.7 


61 


J165112+045917 494.7 


-24.7 ± 4.0 


72 


J172025-005852 1553.3 


37.5 ± 1.3 


234 



Table 7— Continued 



Source P RM SNR 

(mJy) (rad m~^) 



J 172034-005843'' 737.8 49.1 ± 3.7 82 

J184226+794517 316.5 -7.3 ± 4.9 61 

J184150+794728^326.0 -12.0 ± 5.0 60 

J192451-291431 175.5 207.8 ± 11.7 25 

49.1 -738.7 ± 39.2 7.0 

J194114-152431 281.1 -87.2 ± 4.0 74 

J201713+334546 285.3 -388.3 ± 2.5 122 

J211636-205551 220.7 -60.6 ± 10.1 30 

J212344+250410 593.9 -73.1 ± 1.9 152 

J212345+250448 453.4 -72.0 ± 2.6 114 

J222547-045701 290.4 -59.4 ± 8.5 38 

J225357+160853 1002.3 -55.7 ± 1.6 193 

J231956-272713 173.5 -1.0 ± 16.1 20 



^OfF-axis source. 

''Component detected in low resolution 
RM spectrum. 



- 73 - 



With high confidence (7a ~ 40 mJy), we find multiple RM components toward 12 out of 42 
sources, for a total of 61 components. For a lower significance threshold (5cr ~ 30 mJy), we find 
multiple RM components in 19 out of 42 sources, for a total of 92 components. The brightest 
component of each multi-component source ranges from roughly 30% to 90% of the total polarized 
flux in all components, with only a weak dependence on the significance threshold. 

As expected for our frequency coverage, most RM components are unresolved. The multi- 
component RM spectra have components separated by the RM beam size. These multi-component 
spectra are likely caused by how our observing configuration resolves flux from Faraday-thick 
sources. In general, the low-resolution RM spectra have best-flt peaks consistent with the flux- 
weighted mean RM of the components detected at high resolution. This suggests that the low- and 
high-resolution spectra are seeing the same components. This is consistent with the fact that most 
components are thinner than the 50 rad m~'^, high-resolution beam and the low-resolution spectra 
are sensitive to Faraday thicknesses less than 75 rad m~^. 



5. Discussion 

5.1. Peak RM vs. Mean RM 

This survey provides a sample of low- and high-resolution RM spectra to study the effect of 
resolution on measuring RM. Comparing Tables [6] and [7] shows that the mean, low-resolution RM 
is similar to the flux-weighted mean of the components seen at high resolution. In other words, the 
low-resolution spectra are usually smoothed versions of the high-resolution spectra. However, this 
means that the peak RM for a source is resolution dependent. 

Figure HT] shows a histogram of the shift in peak RM between the high- and low-resolution 
spectra in units of significance. An idealized distribution is shown as a Gaussian containing 42 
sources with width la. The observed distribution of RM values has a wider spread than the 
idealized distribution, which refiects the change in peak RM with resolution. 

While the observed distribution of peak RM change is somewhat irregular, it is similar to a 
Gaussian of width 2 — 3a. For a typical RM difference error of 5 rad m~^, this suggests an extra 
RM uncertainty of 10-15 rad m~^ for RM spectra with resolution less than 600 rad m~^. This 
extra uncertainty represents the precision with which one can specify the peak RM of a source in 
our low resolution spectra. 

Generally speaking, any individual RM component cannot be specified more accurately than 
the RM beam width. This is analogous the uncertainty in measuring any centroid of an unre- 
solved cluster of sources. For some ap plications (e.g., measuring mean line-of-sight properties; 



Gaensler et al.ll2005l : iTaylor et al.ll2009l ). the mean RM is adequate, since it averages over all RM 
components. However, any individual component cannot be specified more precisely than the beam 
size due to the uncertainty in the distribution of components within the RM beam. 
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Fig. 47. — Difference in peak RM measured in liigh- and low-resolution RM spectra in terms of the 
significance of the difference. The error in the difference is the quadrature sum of the individual 
RM measurements. The green line shows the observed RM difference and the blue line shows the 
difference expected under Gaussian statistics. 
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Fig. 48. — Left: All RM components measured by the ATA from 1-2 GHz comp ared to |T09| at 



compc 

Toi l 



1.4 G Hz for the 37 central sources. The solid line shows where the ATA and lT09l RMs are equal. 
T09l sources with multiple ATA RM components are connected with a dotted line. The symbol 
size is proportional to the flux density of the component relative to the brightest component in the 
source. Right: A detail of the plot above for small Faraday depths. 



5.2. Comparison with T09 Results 



Figure [48] compares the RM components observed here with the results of lT09l . Since the 
ATA is able to measure multiple significant components per source, the plot shows the relative flux 
in multi-component sources. We find that the fl ux-weighted mean RM of all high-resolution ATA 
components is similar to th e low-resolution, |T09| RM value. The most signiflcant potential difference 
between the ATA and 1t09| values appears in the double source J212344+250410/J212345+250448, 
as disc ussed in ^5.31 This confirms that sources with polarized fiux densities greater than 200 mJy 
in lT09l do not suffer from nvr ambiguity problems and that their RM values are at least 95% reliable 



(more than 35 out of 37 sources). However, as noted in ^5.11 low RM resolution spectra are most 
appropriate for measuring the mean RM along a line of sight, not the peak RM of a source. 

Figure H91 shows the peak RM measured in high-resolution RM spectra (see Table El) projected 
onto the Galactic coo rdinate system. The RM values show large-sc ale structure similar to that seen 
in other RM surveys (|Simard-Normandinlll98ll : iTavlor et al.ll2009l ). Since this map shows the peak 
RM, it is more affected by intrinsic source structure than in low RM resolution observations. 



5.3. Unresolved Double Sources 



As noted in TablelH two pairs of fields in our survey are unresolved by the ATA: (J165111+045919, 
J165112+045917) and (J212344+250410, J212345+250448). These four fields were each observed 
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Fig. 49. — Galactic distribution of the peak RM measured in the high-resolution spectra (Table [6]) 
shown with a Hammer- Aitoff projection. Positive RM is shown with a circle, negative RM is shown 
with a cross. The symbol width scales linearly with the magnitude of the RM, with a minimum 
width similar to the smallest symbols shown. 
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independently, so they provide a unique test of the quahty of ATA RM synthesis analysis. 



The RM spectra toward J165111+045919 and J165112+045917 are consistent with each other. 
Nominally the peak RM is different, but this is a resolution effect, sim ilar to that discussed in 
§5.11 The peak RM measured in two spectra is roughly —5 rad m~ ^ (ITOOl l. This is consistent with 
the flux-weighted mean RM of the indiv idual sources reported by ITOQI . since they are unresolved 
spatially and spectrally. The sum of the |T09| polarized flux in the two components is about twice 
that detected here, so the ATA is beam depolarized by a factor of 2. 

The RM spectra toward two, unresolved parts of 3C 433 (J212344+250410 and J212345+250448) 
are also consistent with each other, showing a single R M com ponent at roughly —72 rad m~^. This 
value differs significantly from th e RM measured by lT09l toward the two sources individually, 
roughly —125 and —94 rad m~^ ( T09l ) . The value seen by the ATA is similar to the va lue of 
—73 lb 1 rad m^^ seen in single-dish, wide-bandwidth observations ( Simard-Normandinlll98lh . The 
low-resolution RM spectra presented here have sensitivity to similar Faraday thicknesses as It09 . 
so their diffe rence (and the consistency of the ATA and single-dish observations) may indicate an 



error m 



T09l . Alternatively, it has been suggested that spatially unresolved pairs of sources can 



confuse the RM synthesis technique in some cases (L. Rudnick, private communication). We cannot 
distinguish between these possibilities with our current analysis technique. 



5.4. Comparison with VLBA RM Images 



Several previous studies have presented milliarc s econd, Very Long Baseline Array (VLBA), po- 
larimetric maps of dozens of AGN (jAsada et al i bnnd : IZavala fc Tavloill2nn3l . [ ' 



2004, 



2005; 



Mantovani et al 



20ld ). These observations were generally done at higher frequencies, ranging from 5 to 22 GHz. 



As the VLBA observations are at higher frequency and have better spatial resolution, they are less 
subject to depth and beam depolarization and show a relatively unambiguous view of the polarized 
emission. Below, we compare the properties observed at high spatial resolution with those of our 
RM spectra. 



5.4.1. 3C 273 



RM maps of 3C 273 (J122906+020305) from 5 to 22 GHz find time-vari able RM with po- 
larized flux density of a few hundred mJy, or a few percent of total intensity (lAsada et al.ll2002l : 



Zavala Sz Tayloiil2005l ). The RM distribution for the core changes from ±3000 rad m ^ over year 



timescales, while the jet (within 10 mas) is predominantly from to 1000 rad m~^. From 1 to 2 
GHz, we find two components at —11 and 151 rad m~^ with polarization fractions of 1 and 0.5%, 
respectively. While 3C 273 has variable RM and brightness distribution, the 1-2 GHz observations 
generally have a lower polarization fraction and RM range than seen b y the VLBA. Another source 



showing this kind of frequency dependence is 3C 446 (see Figure HGJ IZavala &: Tavloijl2003l ). 
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5.4-2. 3C 279 



Our observations of 3C 279 (J125611-054720) find two RM components at 15 and 66 rad m ^ 
with polarized briglitnesses of 270 and 160 mJy (~ 2.7 and 1.6% of total i ntensity). The 1-2 GH z 



RM is similar to the that seen at 8 GHz with the VLBA in June 2000 (IZavala fc Tavloi] 120031) . 
which is much less tha n that seen at earlier dates ( Zavala &: Taylod l200ll ). The observations of 



Zavala &: Taylor! (|2003l ) show two sources at milliarcsecond-resolution with RM of —91 and 13 rad 



m ^ and polarized brightnesses of about 1.4 Jy (~ 9%). Intrinsic variability makes it hard to make 
a strong conclusion, but this similarity suggests that a similar part of 3C 279 is seen from 1 to 8 
GHz. 



5.4.3. 3C 454.3 



Source J225357+160853 is 3C 454.3 (a.k.a. 62251+158), w hich has a low polarizat ion fraction 
of ~ 0.7% at 8 GHz, as compared to 9% in our observations (jZavala &: Tavlori l2003l ) . The RM 
distribution is irregular with a value at peak brightness of —263 rad m~^ at 8 GHz; we measure 
a single component of —54 rad m~^ from 1-2 GHz. This source is a blazar and is highly variable 
at radio wavelengths, so it is difficult to make strong conclusions in this co mparison. However , 
monitoring at 43 GHz shows that its polarization fraction rarely reaches 9% (jJorstad et al.ll2010l ). 
If so, this source may have a higher polarization fraction from 1-2 GHz than at higher frequencies, 
which differs from the trend observed toward 3C 273, 3C 446, and 3C 279. 



5.4.4. M87 

Field J123039+121758 is particularly complicated as it is host to the starburst galaxy M87 
and its jet (see Figur es [22] and [23]) . Comparing the ATA observations to previous work (e.g.. 



Rottmann et al.l Il996l ) shows a similar bright core and faint jet in total intensity. In polarized 
intensity the situation is reversed and the brightest source the ATA sees from 1-2 GHz is the 
termination of the southwest jet; in fact, this polarized source is what was identified by our survey 
selection criteria. In general, the core and other parts of the jet have varying levels of polarized 
flux that change rapidly with frequency. 

The RM spectra of the core and jet of M87 are amongst the most complicated in this survey. 
This is probably not a coincidence, since M87 is one of the few sources that has a core-jet structure 
resolved by the ATA. The spatial separation allows us to resolve multiple Faraday components, 
which otherwise cause depth depolarization. Still, the dramatic increase in polarized flux in the 
1.0 GHz band (especially in the core) suggests that there are strong spectral index effects that 
may confuse the RM cleaning process. Also, the RM values measured for the M87 core in the 
low-resolution RM spectra are very different from those at high-resolution, which is likely related 
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to strong Faraday thickness effects. In principle an algorithm could be developed to disentangle 
these effects, similar to the "multi-frequency synthesis" technique used in synthesis imaging. 

For the M87 jet, we find RM up to ~ 300 rad m~^, while in the M87 core the RM ranges up 
to 100 rad m~^. Despite the relatively large RM we find from 1-2 GHz, M87 h as an even larger 
range of known RM values measured by VLBA imaging {\RM\ < 1000 rad m'^; IZavala fc Tavloil 



20031 ). 



5.5. Limits on Large RM Values 

Limits on RM up to ±9 x 10^ rad m~^ are typically 20 mJy. At the maximum RM, this is 
equivalent to an intrinsic fractional polarization of order 0.5% for our sample. Since this is more 
than an order of magnitude less than the polarization fraction of typical synchrotron sources, these 
components are not likely to be Faraday-thin synchrotron sources. The low-resolution RM spectra 
have a sensitivity to sources with RM thickness up to 75 rad m~^ . If the Faraday thickness is related 
to a physical thickness. Equation [2] requires a physical thickness of 93(ne/l cm~^)~^(B/l fiG)~^ 
pc. However, it is equally likely that the Faraday thickness of the source is caused by multiple RM 
components that are spatially unresolved. 

In general, the 1-2 GHz RM spectra do not detect the large polarization fractions and RM 
values seen by the VLBA from 5-22 GHz. Since the spectral indices for these sources are relatively 
fiat, they should be bright enough to detect from 1-2 GHz. That requires that either (1) the Faraday 
thickness be larger than 75 rad m~^ (in the low-resolution RM spectra), or (2) the polarized source 
structure change between 1 and 10 GHz. 

High-spatial-resoluti on polarimetry of quasars has fourid a d ecreasing polarization fraction 



toward lower frequencies (lAsada et al.ll2002l : iMantovani et al.ll20ld ). which suggests that Faraday- 



thickness depolarization occurs even at high spatial resolution. If this is the sole cause of frequency- 
dependence of the polarization fraction, the lack of large RM requires that Faraday thickness 
increases with Faraday depth. If so, the typical thickness is larger than 75 rad m~^ for depths 
larger than \RM\ > 1300 rad m~^, the largest, high-significance RM detected here. 

It is also likely that wide-b a.ndwidth RM syn t hesis analysis can be confused by relative spectral 



index changes within a source. I Zavala &: Tavloil (|2003l l have shown that the spectral index of the 



3C 454.3 jet is steeper than that of the core, which suggests that 1-2 GHz observations will tend to 
see more of its fiux. The inner core and jet are hosts to the strongest and most organized magnetic 
fields, where one expects large RM and high fractional polarization. The outer jet is likely to be 
interacting with ambient gas, and is thus less organized and has lower RM. The RM synthesis does 
not formally account for these relative changes, so some care is required when interpreting RM 
results from data with large fractional bandwidths. 
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6. Conclusions 



We have conducted a survey of 37 bright, polarized radio sources from 1-2 GHz with the 
ATA. Applying polarimetric calibration allows us to measure the full Stokes information over this 
frequency range. RM synthesis uses these wide-bandwidth spectra to create high-resolution RM 
spectra, sensitive to RM up to 90000 rad m~^ with a resolution as good as 50 rad m"^. Observations 
of 3C 286 show that these data can measure RM with a precision up to 0.6 rad m~^. 

In the 37 fields observed, 42 sources were brightly polarized enough to perform RM synthesis. 
Twelve of the 42 sources have multiple RM components with greater than 7a significance (~ 40 
mJy), for a total of 61 components. Comparing the RM values measured by the ATA at high and 
low RM resolution show that multi-component RM spectra can bias the measurement of the peak 
RM of a source. We show that the peak RM of our sources cannot be known to better than 10-15 
rad m~^ if observed at RM resolution larger than 600 rad m^^. More generally, the peak RM can 
be measured no more accurately than the width of the RM beam. However, for typical applications 
the mean RM measured by low-resolution RM spectra is not affected by this bias. 



The RM spectra measured by the ATA confirm the RM results presented by |T09| for sources 
with polarized flux densities greater than 200 mJy. There is no evidence for rnr ambiguity problems 
and the narrow-bandwidth VLA observations are at least 95% reliable for this sample of sources. 

VLBA maps with milliarcsecond resolution spatially resolve some components that we identify 
in RM spectra. This correspondence suggests that RM synthesis can detect blobs of gas in massive 
AGN jets, which ordinarily must be spatially resolved to be studied. Observations of the time and 
frequency dependence of Faraday thickness and polarization fraction in AGN can be compared to 
simulations to constrain physical models. 

We have shown that RM synthesis is useful because it quantifles the full complexity of the 
Faraday rotation process. Future work will exploit the relative ease of the technique to resolve the 
time-dependent RM structure in AGN. Improving calibration models will also enable wide-field, 
wide-bandwidth RM surveys with the ATA and other radio interferometers. Large, accurate RM 
samples will probe magnetic fields on Galactic and extragalactic scales. 
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